Diffusive instability of a Townsend discharge 

V.V. MikhailenkcE and H.J. Lee 

Pusan National University, 30 Jangjeon-dong, Guemjeong-gu, Pusan 609-735, S. Korea. 

V.S. Mikhailenko 
V.N. Karazm Kharkov National Universtty, 61108 Kharkov, Ukraine. 
(Dated: November 2, 2011) 

The role of the electron diffusion on the stability of a Townsend discharge is investigated. It is 
obtained, that electron diffusion modifies the condition of the steady self-sustenance of the discharge, 
and make discharge unstable. 



I. INTRODUCTION 

There has been considerable interest in non-thermal at- 
mospheric pressure glow discharges over last time due to 
the increased variety of their industrial applications 
The basic feature of the non-thermal discharges is that 
majority of the energy of the applied electric field goes 
into electrons, instead of heating the entire gas in the dis- 
charge cell. Numerous experiments show that depending 
on the parameters of the discharge, atmospheric pressure 
glow discharge is realized in two forms: a Townsend and 
glow discharges. A Townsend discharge is the simplest 
type of glow discharge. It is characterized by the absence 
of quasi-neutral plasma - the absolute value of the ion 
density exceeds much that of the electron density. The 
applied electric field is weakly disturbed by spatial charge 
and the discharge current is governed mainly by the pro- 
cesses of the electron emission from the cathode. The 
current of a Townsend discharge is only limited by the 
external circuit and when the space charge in a Townsend 
discharge becomes large enough to cause a significant dis- 
turbance of the applied field, the transition to glow dis- 
charge occurs. 

The physics of space charge driven transitions of a 
Townsend discharge to subnormal, normal and further 
to abnormal glow has drown considerable attention of 
the discharge investigating community (Refs.[3|-[i3 and 
references therein). Numerous experimental, analytical 
and numerical investigations provide deep insight into 
amazing variety of spatio-temporal processes, which are 
responsible for such transitions. In this paper, the sta- 
bility of a Townsend discharge is investigated for the 
unexplored yet regimes far from such transitions, when 
space charge is too small to produce any significant dis- 
tortion of the applied electric field. We find, however 
that even in that case narrow planar Townsend discharge, 
with the distance between the electrodes considerably 
smaller than the radius of the discharge cell, appears 
unstable. The discovered instability obtains analytical 
confirmation as a resulted from electron diffusion in the 
axial direction. The derivation of the basic equations 



and boundary conditions for the dimensionless variables 
corresponding to the regime of the Townsend discharge 
is presented in Section 2. In Section 3, we present the 
results of the calculations of the modified by electron 
diffusion condition for the steady self-sustenance of the 
Townsend discharge and investigate its stability under 
that condition. A summary and discussion is presented 
in Sec. 4. 



II. BASIC EQUATIONS AND BOUNDARY 
CONDITIONS 

The simplest set of equations containing the basic 
physics necessary for the investigation of the glow dis- 
charge stability comprises the well-known continuity 
equations in the drift-diffusion approximation for elec- 
trons iVe and positive ions Ni, coupled with Poison equa- 
tion for the electrostatic potential $ (e.g.,f6]). 
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Here De, ^J.i, and /Xe are the electron diffusion coefficient 
and mobilities of ions and electrons, respectively; a (E) 
is Townsend's ionization coefficient (e.g., and E is 
the electric field; Eq is the permittivity of free space, 
e is the elementary charge, the Z-axis is directed from 
the cathode to the anode, and t is time. For consid- 
ered here Townsend limit, both mobility and electron 
diffusion coefficient, which otherwise depend on the lo- 
cal electric field, can be assumed as constant. Because 
in the Townsend mode electric field is practically con- 
stant, local field approximation is applicable to ionization 
rate a. Farther, the process of dissociative recombination 
and ion diffusion (we assume that electron temperature 
exceeds greatly the ion temperature) are completely ne- 
glected. In contrast, the axial diffusion of electrons is 
taken into account. The boundary conditions are taken 
in the Townsend form. At the cathode {Z = 0) 

-D,^ + N,^i,E = JN,^i,E, (4) 
oZ 
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where 7 is the secondary electron emission coefficient and 
electrons flax from cathode includes mobility and diffu- 
sive flaxes. At the anode {Z ~ L): 
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and the electric current density at the anode is given by 
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The boundary conditions at the wall of the discharge 
vessel are not relevant for the present one-dimensional 
study. 

In what follows we are interested in processes that are 
assumed to be on time scale which is longer than the ion 
travel time 
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and therefore ion density time derivative will be consid- 
ered as a small perturbation. On that time scale the 
electron density time dependence is eliminated adiabat- 
ically and electron density time derivative in Eg.© will 
be setting to zero. It is convenient to introduce the di- 
mensionless times r, length z and electric field i?, 
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and parameter Se 
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which determines the relative value of the diffusive and 
drift terms in Eq.([2]). By assuming that dimensionless 
anode current is equal to unity, we obtain from boundary 
condition ([9]) the following dimensionless electron density 
rie: 



Similar relation, 
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we use for the dimensionless ion density n^. In dimen- 
sionless form, system (d])-© is 
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where a = aL. Parameter 5, 
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which, as it follows from Eq. p^ . determines the measure 
of the distortion of the external ambient electric field by 
space charge [nl, as well as parameter e^, are assumed to 
be small for a Townsend discharge. 

In dimensionless variables, boundary conditions (|4])- 
([6]) are at the cathode (z — 0) 
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and at the anode {z = 1) 



n. = 0, 
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It is interesting to note, that derivative dE/dr enters 
only into condition ([TS]) at the anode, therefore the 
time dependence of electric field will be determined with 
boundary condition ([T8)) . It follows from Ea.([T4|. that 



E (z, ^) = -S J rii (C, r) rfC + ^0 + 5£ (t) . (19) 

In Eq. we have accounted for that the distortion £ of 
the applied electric field Eq is resulted from space charge. 
It follows from Eq.([19l), that 



(20) 



E{z^I,t) = Eo + 5£{t) 



With electric field (1^ boundary condition at the anode 
becomes 

-l = -ne(z = l,r)(l + Jf (r)) + — , (21) 

dT 

where E^ = l(-Eo = Et) was used. The basic set of 
equations ([T ^ - p^ as well as boundary conditions does 
not contain a time in an explicit form. Therefore in our 
linear stability analysis of the Townsend discharge we 
consider ion and electron densities in a conventional form 

rie.t (z, r; Se, 5) = n^o.^o {z; e^, 6) + Uei^n (z; ee, 5) e^^, 

(22) 

where neo,io is the equilibrium electron (ion) charge and 
nei,ii are their time-dependent small perturbations. Us- 
ing the expansions 

rieOAo {z, r; Se, 5) = n^^^.^o (z; £e) -I- (Jn^o^.o {z; £e) , (23) 
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in Eq. ((2T|) we obtain 

- 1 + 71^0^ (z = 1; Se) + (5n^o = 1; £e) 
= ^ - nfi^ (z = 1; Ee) e^^ - Jn^ = 1; e^) 

CET 

- 5f (r) (nf^ [z = 1; e^) + n'^S = 1; ^e) e^") (25) 

For stationary electron density n'^Q , Eq. (|25)) gives known 
boundary condition 



'eO 

as well as the equation for 8, 



(26) 



(27) 



which, for initial condition E {t — oo) ~ has solution 

8{T)^\-^n^^^{z = l-e,)e^\ (28) 
It also follows from Eq.(|25|) that 

ni^o^(z=l;ee)=0 (29) 

and 
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It stems from Ea. (l30l) . that expansion (|24)) is convergent 
for not small |A|, for which 



<5< lAI. 



(31) 



The system (fT2|) - (fT4|) with boundary conditions (fT6| . 
P7)) . (P5)) - (|5D|) composes the eigenvalue problem for the 
parameter A for the investigations of a Townsend dis- 
charge stability. 



III. DIFFUSIVE INSTABILITY OF A 
TOWNSEND DISCHARGE 



boundary conditions (see Eqs. pJl) . p7)) and ((26l re- 
spectively) nf^ {z = l,ee) = 0, n[,o' (z = l,ee) = 1, 
dnf^^ {z — l,£e) /dz ~ 0, the expressions for electron, 



rigg' (z, Eg), and ion, n^^g'' (z, e^), steady state densities are 
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The condition n''^ (z = 0, £e = 0) = 7n|g'' (z = 0, £e = 0) 
at cathode provides us with condition of the steady self- 
sustenance of Townsend discharge, extended on the ac- 
counting for the effect of electrons diffusion. 



7— (e 
02 



l) = l + 7 



a a2 
( 



(37) 



which is, in fact, the equation which determines Et, used 
in Eq.(|8l) of our transformations to dimensionless vari- 
ables. Last term in Eq.((37|) is negligibly small and it will 
be omitted in what follows. For Ee — Eq. ([37|) reduces 
to well known relationfl] 7 (e" — l) = 1. 

Using presentation ((22|l with S = 0, we obtain the sys- 
tem of equations for ion and electron densities perturba- 
tions, 



9n^°'(z;£e) (0) 



Xnfl' (z; Se) = -anf^ (z; Se) , (38) 



In this section, we solve system (|12l) - p^ with bound- 
ary conditions (|16p . p7)) in asymptotic limit (5 = and 
£e <C 1, for which E = Eq = I. With nomenclature ([^ 
system of equation p2)) - p4)) for for ion, n-g'' (z;£e), and 

electron, n|,g'(z;£e), densities reduces to the following 
simple system 



dn[°^ (z; £e) ^ (0) 



+ an(°g'(z;£e) =0. (33) 
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The solution of Eas. ([32)) - (l33l) provides us with a 
Townsend discharge solution [10|, extended on the ac- 
counting for the effects of electrons diffusion. With 



- £, 
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= anf^^ (z; £e) . (39) 



The solution of that system with boundary conditions 
(113 is 



ni°) (z;£e) = C (^e--(i-^) - ^e--^^"^)) 



(40) 



(32) and 
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where oi 2 are determined by Eg. (1361) and C is arbitrary 
constant of the integration. By using solutions (1391) and 
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(|40)) in boundary condition (IT6l) we obtains the eigenvalue 
equation for parameter A, 
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Accounting for that ai ^ ^ 3> 1, and ai ^ 02 
d (1 — d£e), Eq. p5|) may be simphfied to 



(43) 



(A - aa) (1 - aaEe) = aj {l ~ e""'-^) 



where we assume that A ^ ai or a2. For Se = Eq. p3)) 
has infinite number of complex roots with negative real 
part, which corresponds to the damping oscillations of 
the ion and electron space charges perturbations. Also 
it has two evident real roots, A = 02 and A = 0. Root 
A = 02 is physically senseless, because for X — a2 ion 
density perturbation nf[^ (z; Se) becomes infinite. A = 
corresponds to stationary state, for which presentation 
becomes senseless. This root, however, due to the 
terms with finite Eg in Eq. (j43p . will obtain relatively small 
non-zero value. Expanding the exponential in Ea. (|43p for 
small A, we obtain simple solution for this root. 
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Because the denominator in Ea. (H^ is positive for the 
conditions of a Townsend discharge, root p3)) corre- 
sponds to the growth rate of the aperiodic instability, 
conditioned by electron diffusion. The growth rate (j44|) is 
the main result of this paper. We identify the discovered 
instability as the diffusive instability, because the growth 
rate (|44l) of this instability is proportional to the electron 
diffusion parameter Se- Because of the assumption 6 = 0, 
the instability discovered is true Townsend discharge in- 
stability, which does not lead to the transition to other 
forms of glow discharge. In dimensional physical units, 
the growth rate is equal approximately to 



Ap/iys 



(45) 



where = Te/nii is the ion sound velocity, i^m is ion- 
neutral collision frequency, and Einstein relation, De = 
TejJLe/e, was used. For a discharge in Helium under the 
pressure of 30 Torr, with L = \ mm distance between 



the electrodes, and electron temperature Tg = 
growth rate (|45|) is equal approximately to Xphys 



1 eV, the 
2-10'^ 



IV. CONCLUSION 

In this work, by solving the eigenvalue problem for 
the parameter A, we have shown, that a Townsend dis- 
charge is unstable due to a joint action of the ioniza- 
tion process and electron diffusion. Discovered instability 
is responsible for the spatially homogeneous exponential 
growth with time the perturbations of the electron and 
ion densities due to the ionization processes and from 
this point of view it may be considered as a kind of the 
ionization instability However it is principally dif- 
ferent from known ionization instability of a Townsend 
discharge!^ 01 . The dispersion equation (H^ is derived 
from the boundary value problem solution for which the 
electron diffusion have to be included in the boundary 
conditions (and which can be neglected for the main 
part of the gap between the electrodes). The growth 
rate ()44p is proportional to the ionisation parameter a 
(a2 — a) and does not depend on the local properties of a 
Townsend discharge with practically homogeneous elec- 
tric field, not disturbed by spatial charge, whereas the 
well known ionization instability of a Townsend discharge 
is developed due to violation of the ionization balance re- 
sulted from the local spontaneous fluctuation of the spa- 
tial charge. Accounting for this and that the growth rate 
(|44l) is proportional to the electron diffusion parameter 
Se, we identify the discovered instability as the diffusive 
instability. The presented theory is valid for sufficiently 
small currents and /or strong applied electric field Et, 
for which condition ([31]) is valid. However, because of 
the gradual growth of the ion and electron densities with 
time and concomitant to this process growth of the dis- 
charge current, that in turn leads to the growth of the 
parameter (5, the effects of the small space charge, i5 ^ 1, 
have to be addressed on following stages of the discharge 
development. Therefore, discovered diffusive instability 
may be considered as a precursor of the ordinary spa- 
tially inhomogeneous ionization instability, as well as an 
unstable background for its development at times ■ 
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